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Abstract. We calculate the spin relaxation rates in a parabolic InSb quantum dots due to the 
spin interaction with acoustical phonons. We considered the deformation potential mechanism 
as the dominant electron-phonon coupling in the Pavlov-Firsov spin-phonon Hamiltonian. By 
studying suitable choices of magnetic field and lateral dot size, we determine regions where 
the spin relaxation rates can be practically suppressed. We analyze the behavior of the spin 
relaxation rates as a function of an external magnetic field and mean quantum dot radius. 
Effects of the spin admixture due to Dresselhaus contribution to spin-orbit interaction are also 
discussed. 



1. Introduction 

The ability to manipulate and control processes that involve transitions between spin states is, at 
the moment, of extreme importance due to the recent applications in polarized spin electronics 
and quantum computation. Spin dephasing is the most critical aspect that should be considered 
in the elaboration of proposals of quantum computation based in single spin states as qubits in 
quantum dots (QDs) [1J. While for bulk and for 2D systems the spin relaxation processes has 
been studied in some detail, the problem for QD's still require deeper and further discussions. 
Several processes that can induce spin relaxation in semiconductors have been identified and 
were studied. At the moment remains in discussion which, between these processes, is dominant 
in zero-dimensional systems. Some experimental results have shown good agreement with the 
theoretical predictions for 2D systems [2] but, in general, the identification of the processes 
through direct comparison with the experimental results may become a formidable task. This 
problem is more critical for QDs, since few experimental results exist and the theoretical 
discussion of the spin relaxation mechanisms is still an open subject. Extensive theoretical 
works in QD systems have studied the main phonon mediated spin-flip mechanisms, including 
admixture processes due to spin-orbit coupling [3] and phonon coupling due to interface motion 
(ripple mechanism) [3J. Spin relaxation rates strongly dependent on the dot size, magnetic field 
strength, and temperature, as reported by several authors [3j[5]. It was shown that the quantum 
confinement produces, in general, a strong reduction of the QD relaxation rates. 

In this work, we calculate the spin-flip transition rates, considering the phonon modulation 
by the spin-orbit interaction. For this purpose will use the spin-phonon interaction Hamiltonian 
proposed by Pavlov and Firsov [HE]. In this model, the Hamiltonian describing the transitions 



with spin reversal, due to the scattering of electrons by phonons, can be written in a general 
form, H p h = V p h + (3[a x VVph] • (p + e/cA), where V^h is the phonon operator, a is the spin 
operator, p is the linear momentum operator and A is the vectorial potential related with the 
external magnetic field B. This model has the advantage of being easily adapted to the study 
of other interaction mechanisms with phonons. 



2. Theory 

Based on the effective mass theory applied to the problem of the interaction of an electron with 
lattice vibrations, including the spin-orbit interaction and in presence of an external magnetic 
field, Pavlov and Firsov [6j [7] have obtained the spin-phonon Hamiltonian that describes the 
transitions with spin reversal of the conduction band electrons due to scattering with longitudinal 
lattice vibrations as 
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where, 6 q (6q) are annihilation (creation) phonon operators, the magnetic vector potential A is 
obtained in the symmetric gauge considering an external magnetic field B oriented along the z 
axis, h^ 1 = x ± y, where x, y are unitary vectors along the x and y axis. e q is a unit vectors 
in the direction of the phonon polarization, q is the phonon wave vector, p is the momentum 
operator, v is the average sound velocity, pu is the mass density, V is the system volume and 
d{q) is a coupling constant that depends on the electron-phonon coupling mechanism. Detailed 
expressions for the parameter d(q) can be found in Ref. [7j. 

It has been assumed that the confinement along the z axis is much stronger than the lateral 
confinement. Thus, the lateral motion is decoupled from the one along z and the envelope 
functions separate ip(r) = f(x,y)(j>(z). The z-dependent part of ip( r ) is an eigenfunction of a 
symmetric quantum well of width L. In lens-shaped quasi-two dimensional self assembled QDs, 
the bound states of both electrons and valence-band holes can be understood by assuming a 
lateral spatial confinement modeled by a parabolic potential with rotational symmetry in the 
x — y plane [8], V(p) = ^mujQp 2 , where Uujq is the characteristic confinement energy, and p is 
the radial coordinate. By using the one-band effective mass approximation and considering an 
external magnetic field B applied normal to plane of the QD, the electron lateral wave function 
can be written as 
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where C n j = y/n\/[ir(n + L n is the Laguerre polynomial, n (I) is the principal (azimuthal) 
quantum number, and x{ a ) is the spin wave function for the spin variable a. The corresponding 
eigenenergies are E n ^ a = (2n + \l\ + V)hQ, + (l/2)Tiuj c + (a/2)gp B B, where Q, = (u>q + u^/4) 1 / 2 , 
Hb is the Bohr magneton, a = (Ti/m£l) 1 / 2 is the effective length and uj c = eB/m. In our model, 
we also consider the effects of the Dresselhaus contribution that provides additional admixture 
between spin states. For 2D systems, the linear Dresselhaus Hamiltonian can be written as 
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where p, = —ifiS/i + (e/c)Ai and (3 is the Dresselhaus coupling parameter for this confinement. 
If the confinement potential in the z-direction is considered highly symmetrical, then W z ~ 
and the Rashba contribution can be safely ignored. 

The spin relaxation rates (W) between the electronic states: (n, 1,1 (j)) - * (I)), 
with emission of one acoustic phonon, are calculated from the Fermi golden rule. In the 



Hamiltonian ([T]), we only consider the deformation potential (DP) electron-phonon coupling, 
this is due to the large ^-factor in narrow gap InSb (\g\ ~ 51), the dominant electron-phonon 
coupling for spin relaxation is the DP mechanism [9]. The piezoelectric (PE) coupling governs 
the spin relaxation processes in wide or intermediate gap semiconductors. In the transition 
matrix elements calculation, we not only consider the linear term iq • r in the expansion of 
exp(zq • r) [3], but the integral representation of Bessel function is used in the evaluation of 
electron-phonon overlap integrals. The linear approximation of exp(zq • r) may be valid for spin 
inversion transitions in the spin polarized ground-states of GaAs based QDs where, due to the 
small value of the electron g-factor, only long wavelength phonons are involved. 

3. Results and discussion 

The calculations were performed for a parabolic InSb QD at T ~ K. The material parameters 
for the InSb system are listed in Ref. [10]. We only have considered electron transitions between 
ground state electron Zeeman levels (0, 0, |) — » (0, 0, [) and (0, 1, j) — > (0, 1, |). The temperature 
dependence for one-phonon emission rate is determined from W = Wq(tib + 1), where n# is the 
Bose-Einstein distribution function and Wq is the rate at T = K. In the temperature regime 
T <10 K, we obtain ub + 1 ~ 1 and W ~ Wq. For temperatures larger than few Kelvin degrees, 
two-phonon processes should be considered as the dominant spin relaxation mechanism. These 
types of processes have not been considered in the present calculation. 

In the Figs. [T£t) and [TJa) we show the spin relaxation rates due to DP electron-phonon 
mechanism, as a function of the external magnetic field B and considering some typical values 
for the effective lateral QD size, tq = y/h/mujQ. Some interesting facts about these results should 
be pointed out: i) The rates show a strong dependence with the magnetic field. This fact can be 
explained from the dependence of the rates with the transition energy AE. In general, we obtain 
that W ~ [g* hbB} 11 = (AE) n , n being an integer number that depends on the electron-phonon 
coupling process and g* the effective g-factor. As can be seen in Fig. d]a), when the magnetic 
field increases, the rates also increase until reaching a maximum near B ~ 0.5 T. The position 
of this maximum it is defined from the transition energy conservation: E n i a t — E n ni a i = Uvq. 
ii) The oscillatory behavior of the rates, observed for B > 0.7T are mainly produced by the 
Dresselhaus spin admixture, which modifies the effective Lande <7*-factor. As is showed in Fig. 
Q]a), the g*-factor effects are particulary important for the ground-state Zeeman transition. For 
small magnetic fields, g* — > gbulk an d we may neglect the spin admixture effects. Therefore, the 
spin relaxation shows no oscillations and becomes almost independent of ro- This small QD size 
dependence is in agreement with the experimental observations of Gupta and Kikkawa |llj . iii) 
The rates dependence with the lateral QD size ro, are related to the interplay effects between the 
spatial and magnetic confinements. This competing effects are contained in the electron-phonon 
overlap integral, I oc J f*, v a /(p) exp(iq-r)/ n / i0 -(p)dr. For large fields, the magnetic confinement 
causes a gradual decrease in the overlap integral as the ro increases. For small magnetic fields, 
the spatial confinement is dominant. Thus, when ro diminishes the wave functions become more 
localized and the overlap integral should increase. This effects explain the behavior of the spin 
transition (0, 1, J.) — > (0,1, t) showed in Fig. []] a) (red lines). The Zeeman ground-state rates 
(black lines) are strongly dependent on AE and, for small B, the rates are weakly dependent 
on /. iv) The same rates calculated for GaAs (not showed here), are in general, one order of 
magnitude smaller than InSb rates. As we expected, the relaxation via PE coupling is more 
efficient than via the DP phonon processes. 

In Fig. [T]b) we have plotted the spin relaxation rates for the ground-state Zeeman transition 
as a function of ro and B. We clearly identify a region of strong spin coherence, defined by 
B > 1 T and ro > 100 A. In this regime, the relaxation times are in the ns order and this is an 
important feature for spin qubit engineering. In the B < 0.1T regime, the relaxation times are 
approximately of few /is. This spin frozen region are not robust against the temperature and 
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Figure 1. Spin relaxation rates, W, for a parabolic InSb QD considering the DP coupling 
mechanism. Panel a) shows W as a function of the magnetic field B, for two different electronic 
transitions and several lateral dot radius tq = 50, 75, 100, 125, 150, 175, and 200 A(same tq 
ordering for both transitions), b) Contour plot of the spin relaxation rate as a function of B 
and tq. 



will disappear whenever the thermal energy is larger than the spin transition energy. 
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